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ABSTRACT 

Globular clusters (GCs) are expected to be breeding grounds for the formation of single 
or binary intermediate-mass black holes (IMBHs) of > 100 M Q , but a clear signature 
of their existence is still missing. In this context, we study the process of dynamical 
capture of a millisecond pulsar (MSP) by a single or binary IMBH, simulating various 
types of single-binary and binary-binary encounters. It is found that [IMBH, MSP] 
binaries form over cosmic time in a cluster, at rates < 10 -11 yr , via encounters of 
wide-orbit binary MSPs off the single IMBH, and at a lower pace, via interactions of 
(binary or single) MSPs with the IMBH orbited by a typical cluster star. The formation 
of an [IMBH, MSP] system is strongly inhibited if the IMBH is orbited by a stellar mass 
black hole (BH): in this case, the only viable path is through the formation of a rare 
stable hierarchical triplet with the MSP orbiting exterior to the [IMBH,BH] binary. 
The [IMBH, MSP] binaries that form are relatively short-lived, < 10 8-9 yr, since their 
orbits decay via emission of gravitational waves. The detection of an [IMBH, MSP] 
system has a low probability of occurrence, when inferred from the current sample of 
MSPs in GCs. If next generation radio telescopes, like SKA, will detect an order of 
magnitude larger population of MSP in GCs, at least one [IMBH, MSP] is expected. 
Therefore, a complete search for low-luminosity MSPs in the GCs of the Milky Way 
with SKA will have the potential of testing the hypothesis that IMBHs of order 100 M Q 
are commonly hosted in GCs. The discovery will unambiguously prove that black holes 
exist in the still uncharted interval of masses around > 100 M . 

Key words: Black hole: physics - Globular clusters: general - Stellar dynamics - 
Stars:neutron - Pulsars: general 



1 INTRODUCTION 

1.1 IMBHs: Observations 

A number of observations suggest that intermediate-mass 
black holes (IMBHs) may exist with masses between w 
100 Mq to 10 4 Mq . Heavier than the stellar-mass black holes 
(BHs) born in core-collapse supernovae (3 Mg — 30 Mq; 
Orosz 2003), IMBHs are expected to form in dense, rich stel- 
lar systems through complex dynamical processes. Globular 
clusters (GCs), among the densest stellar systems known in 
galaxies, have therefore become prime sites for their search. 

Gebhardt, Rich & Ho (2002, 2005) suggested the pres- 
ence of an IMBH of 2±o.i x 10 4 M Q , in the cluster Gl of 
M31, on the basis of a joined analysis of photometric and 
spectroscopic measurements. Remarkably, the IMBH in Gl 
seems to lie just on the low-end of the BH mass versus 
one-dimensional dispersion velocity correlation observed in 



spheroids and bulges of nearby galaxies (Ferrarese & Merritt 
2000; Gebhardt et al. 2000). 

In the galactic GC M15, HST and ground-based ob- 
servations of line-of-sight velocities and proper motions, in- 
dicated the occurrence of a central concentration of non- 
luminous matter of 500 j^g Mq , that could be ascribed to 
the presence of an IMBH (van den Bosch et al. 2006; Gerssen 
et al. 2002). By mapping the velocity field, van den Bosch 
et al. (2006) found also evidence of ordered rotation in the 
central 4 arc sec of M15. This unexpected dynamical state 
in a region of rapid relaxation (10 7 yr) may give first evi- 
dence, albeit indirect, that a source of angular momentum 
in the form of a "binary" IMBH may exist in M15 (Mapelli 
et al. 2005). Claims of the possible presence of an IMBH 
have been advanced also in 47 Tucanae (McLaughlin et al. 
2006). 

An additional puzzling picture has emerged from ob- 
servations in the GC NGC 6752. Two millisecond pulsars 
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(MSPs hereon), PSR-B and PSR-E, show unusual accelera- 
tions (D'Amico et al. 2002), that, once ascribed to the over- 
all effect of the cluster potential well, indicate the presence 
of > 1000 Mq of under-luminous matter enclosed within 
the central 0.08 pc (Fcrraro ct al. 2003a). NGC 6752 is 
even more peculiar than M15, since it also hosts two MSPs 
with unusual locations. PSR-A, a binary pulsar with a white 
dwarf (WD) companion (D'Amico et al. 2002; Ferraro et al. 
2003b; Bassa et al. 2003) and a very low orbital eccentricity 
(~ 10~ 5 , D'Amico et al. 2002) holds the record of being the 
farthest MSP ever observed in a GC, at a distance of £s 3.3 
half mass radii. PSR-C, an isolated MSP, ranks second in 
the list of the most offset pulsars known, at a distance of 1.4 
half mass radii from the gravitational center of the cluster 
(D'Amico et al. 2002; Corongiu et al. 2006). Colpi, Possenti 
& Gualandris (2002) first conjectured that PSR-A was pro- 
pelled into the halo in a flyby off a binary BH in the mass 
range between 10 Mq and 100 Mq opening the perspective 
of unveiling binary BHs in GCs (see Section 1.2). Prompted 
by the evidence of under-luminous matter in the core of 
NGC 6752, Colpi, Mapelli & Possenti (2003) carried on an 
extensive analysis of binary-binary encounters with IMBHs, 
to asses the viability of this scenario. They found that a 
~ 100 Mq IMBH with a stellar-mass BH in a binary would 
be the best target for imprinting the necessary thrust to 
PSR-A 1 and at the same time for preserving the low eccen- 
tricity of the binary pulsar (within a factor of 3 for the bulk 
of the simulated encounters). Instead, larger mass IMBHs 
(~ 500 Mq) with star companions can produce the correct 
ejection velocity, but cause the eccentricity to grow much 
larger. Thus, PSR-A had to interact with the very massive 
IMBH only before its recycling phase. 

The observation of IMBHs in GCs is still far from being 
conclusive, since numerical studies have shown that kine- 
matic features as those observed in Gl and M15 can be 
reproduced assuming, in the cluster center, the presence of 
a collection of low-mass compact remnants, with no need of 
a single massive IMBH (Baumgardt et al. 2003a, b) . In addi- 
tion, a single massive (> 1000 M ) IMBH, if present, would 
affect the stellar dynamics (because of energy generation in 
the IMBH cusp) creating a constant density profile of bright 
stars in projection that differs from the typical profile of a 
core-collapse cluster such as M15 (Baumgardt, Makino & 
Ebisuzaki 2004). 

1.2 IMBHs: Theory 

On theoretical ground a number of authors suggested that 
IMBHs may form inside either (i) young and dense star clus- 
ters vulnerable to unstable mass segregation and core col- 
lapse before the most massive stars explode as supernovae 
(Portegies Zwart & McMillan 2002; Freitag, Gurkan & Ra- 
sio 2006; Gurkan et al. 2006) or (ii) dynamically in already 
evolved GCs when all the massive stars have turned into 
stellar-mass BHs (Miller & Hamilton 2002) . In the first case, 

1 Ejection of PSR-A from the core to the halo following exchange 
interactions off normal binary stars can not be excluded, but as 
pointed out by Colpi et al. (2002; Sigurdsson 2003), the binary 
parameters of PSR-A and its evolution make this possibility re- 
mote, and call for fine tuning conditions. 



runaway collisions among young massive stars may lead to 
the formation of a very massive stellar object which ulti- 
mately collapses into an IMBH 2 . In the second case, IMBH 
formation requires a succession of close gravitational en- 
counters among stellar-mass BHs: being the heaviest objects 
in the cluster, these BHs may segregate in the core under the 
action of the Spitzer's mass stratification instability (Spitzer 
1969; Lightman & Fall 1978; Watters, Joshi & Rasio 2000; 
Khalisi, Amaro-Seoane & Spurzem 2006), forming a dense 
core which becomes dynamically decoupled from the rest of 
the stars. Hardening and recoil among the interacting BHs 
lead to their ejection from the cluster (Sigurdsson & Hern- 
quist 1993; Kulkarni, Hut & McMillan 1993; Portegies Zwart 
& McMillan 2000) and at the same time to the increase of 
their mass because of repeated mergers (Miller & Hamil- 
ton 2002). O'Leary et al. (2006) have recently shown that 
there is a significant probability (between 20% to 80%) of 
BH growth, and found final masses > 100 Mq. After evap- 
oration of most of the BHs on a timescale of ~ Gyr, one 
IMBH and/or few BHs, single or in binaries, may remain 
inside the GC. 

The recent discovery of a luminous, highly variable X- 
ray source in one GC of NGC 4472 (Maccarone et al. 2007) 
may have just provided first evidence that at least one BH 
is retained inside. Whether this source in NGC 4472 is an 
accreting BH or IMBH is still uncertain, but this finding 
goes in the direction noted by Pfahl (2005), who considered 
the possibility that an IMBH would tidally capture a star 
leading to the turn-on of a bright X-ray source. 

Given all these uncertainties and the importance of es- 
tablishing the possible existence of IMBH in GCs, we explore 
in this paper an alternative root, i.e., the possibility that 
gravitational encounters off the IMBH provide a path for 
the dynamical capture of a MSP and the formation of a bi- 
nary (hereafter labeled [IMBH,MSP]) comprising the IMBH 
and the MSP. Timing of the radio signal emitted by the MSP 
would provide in this way a direct, unambiguous measure of 
the BH mass. 

Motivated by the observation of the halo MSPs in NGC 
6752, we simulate a series of dynamical interactions between 
a binary MSP and a single or a binary IMBH, and also 
between a single MSP and a binary IMBH. In the context 
adopted, the binary IMBH may have a stellar-mass BH, or 
a star, as companion. 

The outline of the paper is as follows. In Section 2, we 
describe the initial conditions of the three and four-body 
encounters. In Section 3, we compute cross sections for the 
formation of [IMBH, MSP] systems coming from encounters 
with PSR-A like MSP binaries. We study the orbital charac- 
teristics of the [IMBH, MSP] binaries in their end-states, and 
explore the stability of triple systems that form, against dy- 
namical and resonant self-interactions. Binary systems com- 
posed by the WD and the IMBH are also considered, and 

2 The effects of the environment, of rotation and metallicity, on 
the formation and fate of these ultra— massive stars are largely 
unknown. A recent study on the mass loss of merged stars (during 
and after the merger) of ~ 100 Mq have shown that this does not 
seem to inhibit the formation of very massive stars (Suzuki et 
al. 2007). However further studies are needed in order to better 
constrain the evolution of those more massive object (~ 1000 Mq) 
that should form ~ 1000 M IMBH. 
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Table 1. Initial parameters for simulations with PSR-A like ini- 
tial MSP binaries. Rows refer to different initial states of the 
IMBH (referred as channels in the text). The different columns 
refer to: selected IMBH mass, minimum and maximum values for 
the distribution of the semi-major axis (for the [IMBH, star] and 
[IMBH,BH] binaries) and number of runs for each simulation. The 
first eight lines refer to encounters with the [MSP, WD] binary, the 
last two refer to encounters with a single MSP. 



the results are shortly summarized in Section 4. In Section 
5, we show the results obtained from simulations with bi- 
nary MSPs different from PSR-A that represent the ob- 
served population in GCs. We study their end-states and 
their characteristic lifetimes taking into account for their 
hardening by cluster stars and by gravitational wave driven 
in-spiral. In Section 6 we study the detectability of MSPs 
around IMBHs in GCs and discuss the potential importance 
of these systems for next-generation deep radio surveys in 
the Galactic halo. In Section 7 we summarize our findings. 



2 GRAVITATIONAL ENCOUNTERS 
2.1 The projectile 

We consider encounters in which the projectile is either a 
[MSP, WD] binary, or a single MSP. As first case-study, 
we simulate [MSP, WD] systems similar to PSR-A in NGC 
6752: the MSP has a mass ttimsp = 1.4 M Q and a WD com- 
panion of mwD = 0.2 M Q ; the binary has semi-major axis 
«MSP,i = 0.0223 AU, orbital period of 0.86 days, and orbital 
eccentricity eMSP.i = 10~ 5 . 

We then simulate binary MSPs whose characteristics 
are extrapolated from the observed sample of MSPs belong- 
ing to the GCs of the Milky Way (Camilo & Rasio 2005) 
(see Section 5 for further discussion). For the single MSP, 
we consider ttimsp = 1.4 Mq. 



2.2 The target IMBH 

The target is an IMBH, either single or binary, and has 
no stellar cusp (Baumgardt et al 2004). In agreement with 
O'Leary et al. (2006) and Colpi et al.(2003), its mass Mimbh 
is either 100 M Q or 300 M . 



The binary IMBHs have initial semi-major axes and 
eccentricities drawn from probability distributions that ac- 
count for their physical conditions in a GC. In details, the 
initial properties of the target [IMBH, star] and [IMBH, BH] 
binaries are the following. 

• [IMBH, star]: We randomly generate the mass m» of 
the star, the semi-major axis a* and the eccentricity e* . The 
values for m* follow a current mass function biased toward 
massive stars, in order to account for dynamical mass seg- 
regation in the core of the cluster. We thus consider a mass 
function dN /dm oc m~( 1+x ) with x = —5 as inferred from 
observations of 47 Tucanae (Monkman et al. 2006) with an 
upper cut-off mass of 0.95 M Q . For the semi-major axes 
we follow the analysis proposed by Pfahl (2005) and briefly 
summarized in Appendix A. The values of a* refer to con- 
ditions acquired in dynamical ionization of incoming stellar 
binaries off an initially single IMBH. Table 1 gives the initial 
minimum and maximum semi-major axes used at the onset 
of the simulations. The eccentricity e* follows a thermal dis- 
tribution (Blecha et al. 2006). The same distribution for o*, 
e* and m» is used for the interaction of the [IMBH, star] bi- 
nary both with [MSP, WD] and single MSP. To distinguish 
these two cases, hereon we will refer to the latter using the 
subscript "MSP,single" . 

• [IMBH, BH] : The IMBH has a BH companion of m BH = 
10 Mq . The binary has semi-major axis ap,H drawn from two 
distinct probability distributions, which have been derived: 
(i) from the hardening due to encounters off cluster stars 
(subscript [h,*], hereon), occurring on a time-scale (Quinlan 
1996; Mapelli et al. 2005) 
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where (p*), («*) and £ are the mean stellar mass density, dis- 
persion velocity and hardening efficiency (we assume (v*) = 
lOwiokms -1 , 5 = 1 (Colpi et al. 2003), obh = 5as AU and 
for the density (p*) = 7 x 10 5 pr,.» M© pc~ 3 , the value in- 
ferred averaging over the GC sample currently hosting the 
population of known MSPs (see Section 5)); 
(ii) from the in-spiral driven by gravitational wave back- 
reaction (subscript [gw], hereon), when the binary is tight 
(Section 1 of Appendix A, for details). The corresponding 
time-scale, function of the semi-major axis <zbh and eccen- 
tricity eaH (Peters & Mathews 1963), is: 
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where the following normalizations are used to estimate 
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0.7: 



a BH = 0.2o .2 AU, Mimbh = 
100MiooM Q , m BH = lOmi O M , and M T = Mimbh + 
tobe = IIOMt.hoMq. The peak of the composite semi- 
major axis distribution occurs when the two processes be- 
come comparable, i.e. at a distance 
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corresponding to th = i g w, inferred from equations (1) and 
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(2). Typical separations for our [IMBH,BH] binaries are ~ 
0.3 AU. 

In the hardening phase by stars the eccentricity cbh is ex- 
tracted from a thermal distribution, while during the grav- 
itational wave driven phase the values of eBH are inferred 
considering the modifications induced by gravitational wave 
loss (see Section f of Appendix A) . 



2.3 Code and outcomes 

We run the numerical code Chain (kindly suited by S. 
Aarseth) which makes use of a Bulirsch-Stoer variable step 
integrator with KS-chain regularization. The code FEBO 
(FEw-BOdy), based on a fifth-order Runge-Kutta scheme 
(described in Colpi, Mapelli & Possenti 2003 and in Mapelli 
et al. 2005), has been used for trial runs and gives results in 
nice agreement with Chain. 

The impact parameters of the incoming binaries are dis- 
tributed uniformly in b 2 (Hut & Bahcall 1983) up to a max- 
imum value &max ( see Section 3 of Appendix A) . The phases 
of the binaries and the angles describing the initial direc- 
tion and inclination of the encounter are extracted from the 
distributions by Hut & Bahcall (1983). The relative speed 
t>oo has been sampled at random from a uniform distribu- 
tion, in the range 8-12 km s~ , consistent with the values 
of NGC 6752 (Dubath, Meylan & Mayor 1997). The rela- 
tive distance between the centers of mass of the interacting 
binaries is set equal to the gravitational influence radius of 
the target IMBH, r inf ~ 2GMimbh/{i)oo) 2 (~ 2000 AU for 
the 100 Mq case 3 , obtained for a stellar dispersion of 10 km 

After each single-binary encounter we can classify the 
end-states as following: 

(A) Fly-by: the binary maintains its components, but it can 
exit with a different energy and angular momentum; 

(B) Tidal disruption: the interacting binary is broken by the 
massive IMBH. The tidal disruption can end with an ioniza- 
tion (B.l), if the final system consists of three single bodies, 
or with an exchange (B.2), if one of the two components 
is captured by the single. The tidal perturbation occurs at 
a distance r T = a M SP,i [A^imbh /(jtimsp + 7«wd)] 1/3 , where 
the gradient exerted by the IMBH on the incoming binary 
exceeds its binding energy. For our binary pulsar, tt ~ 0.1 
AU. 

In the case of binary-binary encounters the possible end- 
states are analogous (i.e. fly-bies and tidal disruptions), but 
complicated by the fact that the interacting binaries are two. 
In particular, we can observe the tidal disruption of only one 
of the two binaries (mostly the softer [MSP, WD] binary), or 
of both of them. After the tidal disruption of the [MSP, WD] 
binary: 

(B.l) The [MSP, WD] can be fully ionized (i.e. both compo- 
nents escape); 

(B.2) One of the two components remains bound to 



3 For the 300M Q IMBH, the larger initial distance (6000 AU) 
makes prohibitive the integration time for the simulations run 
with FEBO. For this reason integration starts at 2000 AU after 
correcting for the relative parabolic motion. For consistency, wc 
have chosen to adopt the same corrections also for the simulations 
run with Chain. 



the [IMBH, star] or [IMBH,BH] binary, forming a sta- 
ble/unstable triplet. Some triplets show a characteristic con- 
figuration of two nested binaries, where two of the three 
components are bound in a tight binary, while the other one 
orbits around. This type of systems are termed hierarchical 
triplets. 

A hierarchical triple is stable if it satisfies the relation 
(Mardling & Aarseth 1999) 



R, 
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where R p is the pericenter of the outer binary, a ln the semi- 
major axis of the inner binary, e ou the eccentricity of the 
outer binary and q the mass ratio between the external com- 
ponent and the inner binary. If the triplet is unstable, the 
evolution of the system ends with the expulsion of one of 
the three components (preferentially, the less bound com- 
panion) . 

In the simulations, the integration is halted when the 
outgoing unbound star(s) is (are) at a sufficiently large dis- 
tance from the center of mass of the target binary or of the 
newly formed binary (or triplet). This maximum distance 
has been chosen equal to 50 times the semi-major axis of 
the system left. If the outgoing star (or binary) is still at 
such a distance after at least 2000 time-units, we stop the 
integration and we classify the encounter as an unresolved 
resonance. 



3 [IMBH,MSP] BINARIES 
3.1 Cross Sections 

We are interested in deriving the frequency of encounters 
ending with the formation of a [IMBH, MSP] binary. Thus, 
we computed fx = Nx/N, i.e. the probability factor asso- 
ciated to channel X, where A is the total number of runs, 
and Ax is the number of cases in which event X occurs. The 
cross section for channel X can be written as 



S x = tv fx b n 



(5) 



where 6„ al is the square of the maximum impact parameter 
that includes "all" relevant encounters leading to X (Sig- 
urdsson & Phinney 1993; see Section 3 of Appendix A for 
its operative definition). Table 2 summarizes our results. 

• In the encounters between the single IMBH and the 
[MSP,WD], we find that ionization of the incoming binary 
leads to the formation of [IMBH, MSP] systems with an oc- 
currence ~ 10%. The cross section in physical units is about 
a few hundreds AU 2 , and increases with the IMBH mass. 

• [IMBH, star]: In the case of binary-binary encounters 
with the target binary [IMBH, star], we often observe the 
exchange between the star and the heavier MSP, leading 
to the formation of an [IMBH, MSP] binary. The cross sec- 
tion for the formation of the [IMBH, MSP] binary is slightly 
larger than for the isolated IMBH in the case of an IMBH of 
100 M Q , whereas the opposite holds for an IMBH of 100 M . 

• [IMBH, star]MSP single: In the encounter 
of the [IMBH,star] and the single MSP we again observe 
the exchange of the star with the MSP, thus forming an 
[IMBH, MSP] system. We note that the frequency is a factor 
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[100, star] M SP, single 1-6 - 126 



[300, star] M sp, single 0.65 - 66 



Table 2. Occurrence fractions (/msp an d /wd) and cross sections (Smsp an d ^Wd) calculated from equation (5) of [IMBH,MSP] and 
[IMBH, WD] binaries for each initial state of the IMBH, and for PSR-A like MSP binaries. Bracket (tr,in) denotes the occurrence of 
stable triplets where the MSP or the WD binds forming the inner binary. Bracket (tr,ou) denotes the occurrence of stable triplets where 
the MSP or WD binds forming the outer binary. The last two lines correspond to exchanges of a single MSP off the [IMBH, star] binary. 
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Table 3. Probability coefficient tux as defined in Section 6, rates of formation of observable [IMBH, MSP] binaries, and lifetimes tufe MSPi 
for (v*)=10 kms — 1 , (p„) = 7 X 10 5 Mq pc~ 3 . The channels of formation are the same as in Tabic 1. 

somewhat lower for the single MSP than in the [MSP, WD] 
case and this involves smaller cross sections too. 

• [IMBH,BH]: In general, the presence of a massive com- 
panion such as a stellar-mass BH does not favor the for- 
mation of an [IMBH, MSP], since the exchange probability 
is negligible. Triple systems may alternatively form. In rare 
cases (< 0.1%) stable triplets can form with the MSP mem- 
ber of the inner binary [(IMBH,MSP) ,BH] . This occurs when 
the IMBH binary is in its hardening phase by dynamical en- 
counters. When the [IMBH,BH] is in the phase of hardening 
by emission of gravitational waves, the MSP binds to the 
[IMBH,BH] as external companion with an higher probabil- 
ity (fx ~ 0.2 — 0.3%) than in the hardening by scattering 
regime. 

3.2 [IMBH,MSP] binary parameters 

In this section we explore the properties of the [IMBH, MSP] 
systems that have formed dynamically. Fig. 1 shows the dis- 
tribution of semi-major axes resulting from encounters with 
the 100 M Q IMBH. In the case of tidal disruption of the 



[MSP, WD] off the single IMBH, we find that the distribu- 
tion peaks at ~ 1 AU. This value agrees with the analytical 
estimate (Pfahl 2005) obtained in the impulse approxima- 
tion, i.e. considering that the incoming [MSP, WD] binary 
is approaching the IMBH along a parabolic orbit, and that 
is disrupted instantaneously at the tidal radius tt- Accord- 
ing to this analytical model (Pfahl 2005), the most likely 
end-state has a binding energy per unit mass 

E rrvwD ^ 

m M sp + mwD 

where Vt ~ (GMimbh/'*t) 1 ' /2 and Kei is the relative veloc- 
ity of the [MSP, WD] binary before the encounter. The corre- 
sponding semi-major axis of the newly formed [IMBH, MSP] 
binary is 

aMSP.i Mimbh / rriMSP + m W D \ 1/3 ,_, 

OMSP.f ~ - 7= I (0 

2V2 mwD V Mimbh / 

which perfectly agrees with the results of our simulations 
(oMSP.f ~ 1 AU for a 100 M IMBH). 

Fig. 1 also shows the distribution of the semi-major 
axes of the [IMBH,MSP] formed during the [MSP, WD] in- 
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Figure 1. Distribution of the semi- major axes of the 
[IMBH, MSP] binaries, normalized to the corresponding fraction 
of events, for PSR-A like initial MSP binaries. The IMBH has 
a mass of 100 Mq. Shaded histogram with dotted lines refers 
to [IMBH, MSP] systems formed after tidal disruption off the 
single IMBH. Shaded histogram with solid lines refers to the 
[IMBH, MSP] binaries that form after the exchange of the initial 
star in the [IMBH, star] binary. 
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Figure 2. Distribution of eccentricities of [IMBH, MSP] binaries, 
normalized to the corresponding fraction of events. Shaded his- 
tograms refer to the same cases as in Fig. 1. 



teraction off the [IMBH, star] binary, following the disrup- 
tion of the [MSP, WD] at ~ n and the subsequent exchange 
of the MSP off the star. The MSP is captured on a close 
orbit, and, from simple energy arguments, the most likely 
end-state is expected to have a specific energy 



mwD 



mMSP + mwD 



-VtV« 



m, GMimbh 
a» 2mMSP 



(8) 



Indeed, during the triple encounter between the MSP, the 
star and the IMBH (after the expulsion of the WD), an 
energy (at least) equal to the binding energy of the star 
before ejection needs to be extracted, in order to unbind the 
star. The characteristic semi-major axis of the newly formed 
[IMBH,MSP] will thus be 

aMSP.f 



a MSP,f 



1 + (m./miisp) a M SP,f/«* 



(9) 



If we consider mean values for the initial m*/a* selecting 
all the systems that end with an [IMBH, MSP] binary, we 
find m*/a* ~ I.68M0/AU. This corresponds to a semi- 
analytical estimate aMSP.f ~ 0-45 AU, in good agreement 
with the peak of the corresponding semi-major axis distri- 
bution derived from our simulations (Fig. 1). 

Fig. 2 shows the distribution of the eccentricities for the 
same binaries. For the case of tidal capture the eccentrici- 
ties at which the MSP binds to the IMBH are above 0.9; 
for the formation channel through exchange the spread of 
the final eccentricity distribution is much larger, according 
to a thermal distribution. This can eventually be the effect 
of repeated interactions between the MSP and the initial 
companion of the IMBH during the transient state of un- 
stable triplet. The distribution of the semi-major axis and 
eccentricity of [IMBH,MSP]msp single systems formed by the 
exchange off the single MSP are similar to the ones formed 
in the interaction of the [MSP, WD] off the [IMBH,star]. 

Finally we note that in the case of a 300 Mq IMBH, the 
distributions are similar and only slightly skewed to larger 
values of the semi-major axes, as should be expected for a 
more massive BH (see equation 7). 



3.3 Hierarchical triplets 

As previously noted, the only way a MSP can be retained in 
the presence of an [IMBH,BH] binary is through the forma- 
tion of hierarchical stable triple systems. Two possibilities 
exist: either the formation of a [(IMBH,MSP),BH] where the 
MSP is closely bound to the IMBH, or the formation of a 
[(IMBH,BH),MSP] with the MSP as external object. 

Triple systems of the first type are rare, because the 
MSP tends to bind preferentially on orbits where its motion 
is gravitationally perturbed by the stellar-mass BH causing 
the MSP to be finally ejected. Only triplets of the second 
type are seen to form with a non negligible probability (~ 
0.2%): the MSP binds on very wide (20-100 AU), eccentric 
orbits (> 0.6), as shown in Figs. 3 and 4. The triplets in 
consideration are extremely hierarchical (i.e., -Rmsp.ou 3> 
ftBH.in), in order to fulfill the stability condition. 

Hierarchical triplets of this type are likely to survive 
inside the GC and to turn into a [IMBH, MSP]. Indeed, once 
the triplet has formed, the MSP shrinks its orbit with time 
due to dynamical encounters off cluster stars while the inner 
binary hardens due to gravitational wave emission. Since the 
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Figure 3. MSP semi-major axis omsp.ou of the outer binary 
versus semi-major axis bbh in of the inner binary (IMBH,BH) 
of stable hierarchical triple systems. The plot refers to an initial 
[IMBH,BH] binary of 100 M and 10 M , and a initial PSR-A- 
like MSP binary. 



Figure 4. MSP eccentricity eMSP,ou of the outer binary versus 
eccentricity cjh ; n of the inner binary (IMBH,BH) of stable hi- 
erarchical triple systems: the initial parameters of the involved 
binaries are the same as in Fig. 3. 



hardening time of the inner binary is usually shorter than 
that of the outer binary, these triplets are transient states 
ending with the formation of a new [IMBH,MSP] binary 
following BH coalescence. 



4 [IMBH,WD] BINARIES 

For the sake of completeness, the results on the formation 
of [IMBH, WD] binaries are also summarized in Table 2. In 
the case of the capture of the WD by the single IMBH, 
we note that the occurrence fraction of [IMBH, WD] is only 
slightly lower than that of [IMBH, MSP] while it decreases 
of a factor ~ 5 for the [IMBH, star] cases, as shown in Table 
2. If the IMBH has a companion star, the WD preferentially 
binds in triplet configurations. In fact the WD can be re- 
tained around the IMBH only if it forms a hierarchical triplet 
[(IMBH, WD), star]. This is due to the smaller mass of the 
WD relative to the star that makes exchanges very unlikely. 
The same is true for the [IMBH,BH] cases: stable triplets 
form with the WD in the inner binary, i.e [(IMBH,WD) ,BH] , 
when the IMBH binary is hardening by scattering stars. On 
the contrary, the fraction of stable triplets significantly drops 
during the gravitational wave driven phase (~ 0.04%). This 
is due to the fact that the WD preferentially binds to the 
IMBH on a orbit strongly perturbed by the stellar mass BH. 
The cross sections computed using equation (5) are reported 
in Table 2 and their values reflect their dependence upon fx- 
Fig. 5 shows the distributions of the semi-major axis 
and eccentricity for the WD case, considering only the in- 
teraction with the single IMBH. Because of its lighter mass 
with respect to the MSP, the WD binds around the single 



IMBH on tighter orbits and the peak is around 0.17 AU, in 
agreement with Pfahl's analysis (2005) 4 . 

The channel that we have outlined for the formation 
of a [IMBH,WD] binary is probably not the dominant one, 
because of the higher number of [WD, star] with respect to 
[MSP, WD] binaries. For this reason we have chosen not to 
discuss the formation rate of [IMBH, WD] binaries in more 
details. 



5 [IMBH,MSP] IN GLOBULAR CLUSTERS 

So far, we have considered only binary MSPs which mimic 
the properties of PSR-A in NGC 6752. Compared to PSR-A 
however, binary MSPs in GCs display a wider distribution of 
properties in their orbits and masses (Camilo & Rasio 2005). 
Since the cross section for the formation of [IMBH, MSP] 
systems as well as their ending states depend on the initial 
semi-major axes and total mass of the impinging [MSP, WD] 
binaries, in this section we have simulated a set of interac- 
tions varying the properties of the binary MSP. 

Binary MSPs in GCs show a double peaked distribution 
of their semi-major axes in the interval [0.0024 AU, 0.035 
AU], while a number of "outliers" spread over larger orbital 
separations (see Fig. 3 in Camilo & Rasio 2005). Outliers 
count for the 25% of the entire population. We have fitted 
the observed distribution with (i) an asymmetric Landau 

4 If the WD is captured instead of the MSP, equation (7) is mod- 
ified to take into account for the different mass of the expelled 
star, thus giving a WD , ~ ( "C" )'^ = 
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Figure 5. Distributions of semi-major axis and eccentricities of 
the [IMBH, WD] binaries, normalized to the corresponding frac- 
tion of events, for the single IMBH of 100 M Q , and a PSR-A like 
initial MSP binary. 



profile, peaked at 0.005 AU, in the range [0.0024 AU, 0.02 
AU] (defining class I [short period binary MSPs]), plus (ii) 
a Gaussian profile, centered around 0.026 AU, in the range 
[0.02 AU, 0.035 AU] (defining class II [long period binary 
MSPs]). According to Camilo & Rasio (2005), we have as- 
signed a companion WD mass of 0.03 Mq for class I, and 
of mwD = 0.2 Mq for class II. For the binary MSPs re- 
ferred to as outliers, we have taken omspj = 0.21 AU and 
mwD = 0.34 Mq, corresponding to their mean properties. 



5.1 Cross sections 

Table 4 collects the results obtained considering as target an 
IMBH of 100 M . We find, in the case of the single IMBH, 
that the cross section is larger for the outliers compared to 
class I+II, due to their initially wider separation. For the 
[IMBH, MSP] binaries formed following the exchange of the 
initial stellar companion we obtain similar results, but the 



Figure 6. Distribution of the semi-major axes of [IMBH, MSP] 
binaries, normalized to the corresponding fraction of events. The 
IMBH has a mass of 100 Mq. Left panel refers to encounters off 
the single IMBH; solid, dotted and dot-dashed, lines refer to scat- 
tering with class I, class II and outliers, respectively. Right panel 
refers to encounters off the [IMBH, star] binary: solid, and dashed 
lines refer to class I+II, and outliers, respectively. 



differences in cross section between outliers and class I and 
II is less pronounced. 



5.2 Orbital parameters 

Fig. 6 (left panel) shows the distributions of the semi-major 
axes of the [IMBH, MSP] binaries formed after the interac- 
tions off a single IMBH. It appears that different populations 
of [MSP, WD] binaries lead to the formation of [IMBH,MSP] 
systems with different orbital characteristics. The peak of 
the semi-major axis distribution for each class can be in- 
ferred from equation (7): 1.7 AU for the short period, class 
I binaries, 1.1 AU for the long period, class II binaries, and 
5.6 AU for the outliers. A clear trend is also visible for the 
eccentricities (Fig. 7 left panel): the lighter the WD is, the 
more eccentric (and with a narrower spread) is the orbit of 
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Table 4. Outcomes from the encounters of different kinds of binary MSPs in GCs with a single or a binary IMBH of 100 Mq. Columns: 
number N of runs for each set of simulations , occurrence fraction (/msp normalized to N), cross section Smsp ( as defined in Section 
3.1), probability coefficient wx as defined in Section 6, characteristic formation rates Tx, and lifetimes tm e (estimated as in Section 5.3). 
These times are computed considering (v t )=10 kms" 1 , (p„) = 7 X 10 5 Mq pc — 3 , and a core radius of 0.75 pc. First (last) two rows refer 
to encounters with class I+II binaries and to outliers scattering off a single (binary) IMBH, respectively. 
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the [IMBH, MSP] binary. This correlation is due to angular 
momentum conservation: 



mwD 



GaMSP.i 

m M SP + mwD 



M 1MW J GaMSpAl - e <> 



rriMSP + Mimbh 



Using equation (7) this implies 

i 2 3 / \-4/3 

1 — e f oc m WD (m M sp + m WD ) 



(10) 



(11) 



Fig. 6 (right panel) shows the semi-major axes of the 
[IMBH, MSP] systems formed after the interaction with 
the [IMBH, star] systems. The distributions are skewed to 
smaller separations, compared to the case of a single IMBH, 
due to the fact that the MSP has ejected the star (see Sec- 
tion 3.2). The smaller cross section for the [IMBH,star] case 
compared to the single IMBH, for the family of the outliers 
(see Table 4), is due to the occurrence of unstable triplets 
where the MSP, that binds onto wider orbit (see equation 
(7)), is preferentially expelled. Fig. 7 (right panel) shows 
the eccentricity distribution, relative to encounters off the 
[IMBH,star] binaries, which it turns out similar to that of 
Fig. 2. 



5.3 Lifetimes 

The simulations provide the semi-major axes and eccentric- 
ities of the [IMBH, MSP] systems formed. So, using equa- 
tions (A2), (A3) and (A4) of Section 1 of Appendix A, we 
can calculate their subsequent orbital evolution, controlled 
either by hardening off cluster stars or by gravitational wave 
back-reaction. The lifetime is defined as the sum of the time 
necessary for the individual binary to harden by stars until 
the separation a gw (equation (3)) is attained, plus the time 
for gravitational wave in-spiral at a gw , i.e., tnt e = ih + igw 
The mean values of the binary lifetimes are reported in Table 
3 for PSR-A-like initial MSP binaries, and in Table 4 for the 
complete population. Note that tnte is computed assuming 
that the eccentricity cmsp does not vary during the hard- 
ening phase against stars. A further increase in cmsp can 
bring the binary into the gravitational waves regime faster, 
while a reduction can make the binary more long-lived. The 
[MSP,IMBH] binaries formed are already very eccentric. If 
dynamical interaction tends to bring the eccentricity distri- 
bution closer to the thermal one, we then can argue that our 
estimated lifetimes represent lower limits. 

Fig. 8 shows the characteristic lifetimes of the 
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[IMBH, MSP] binaries described in Section 5.2. Left panel 
refers to encounters off the single 100 M IMBH. We note 
that the different families of [MSP, WD] binaries lead to 
[IMBH, MSP] systems with different lifetimes: in particular 
for class I, tn{ e ~ 6 x 10 7 yr due to the extremely high 
eccentricities at which the new systems form. By contrast, 
class II and the outliers have higher t\it e > 10 8 yr. Right 
panel of Fig. 8 refers to encounters of MSP binaries off the 
[IMBH,star] system. In this case the distributions seem not 
to depend strongly on the incoming binaries: outliers as well 
as class I and II show very similar lifetime distributions with 
characteristic values around 4 x 10 8 yr. 



and estimate the number of expected systems in the Milky 
Way GCs. 

The rate of formation for channel X reads 



Tx ~ JlMSPWx(«oo}Sx 



(12) 



where timsp is the number density of MSPs in the cluster 
core of radius r c , Ex the cross section defined in equation 
(5) and tux the probability coefficient (estimated below), 
associated to channel X. 

The structural parameters of GCs span a large interval 
of values. In order to estimate Tx, we considered only the 
23 GCs that are known to host at least one MSP. For each 
GC in this selected sample, we computed the MSP number 
density as timsp ~ A r MSp/4rj? where A^msp is half of the 
number of currently observed MSPs in every GC in order to 
take into account the fact that not all MSPs are hosted inside 
the GCs core. The mean value of umsp obtained considering 
the sample of galactic GCs is « 2 x 10~ 14 A IP 3 . 

For the calculation of wx, we adopted a ratio of 2 for the 
relative number of single and binary MSPs, in accordance 
with the ratio observed (Camilo & Rasio 2005). The outliers 
account for 25% of the binary MSPs, and class I and II for 
50% and 25%, respectively. Following Blecha et al. (2006), 
we also assume that the IMBH lives as single object for 
~ 40% of its lifetime, whereas for the remaining ~ 60% it is 
bound with a cluster star. The values of iox are computed 
according to these simple recipes and are collected in Tables 
3 and 4 together with the estimates mean rates Tx- We 
note that the main contribution comes from binary MSPs 
belonging to the family of the outliers, scattering off the 
single IMBH. 

As previously discussed in Section 5.3 and shown in 
Fig. 8, the [IMBH, MSP] binaries have characteristic life- 
times shorter than their typical formation timescales. Conse- 
quently, the expected number of [IMBH, MSP] binaries that 
formed and reside in a GC is roughly given by 



Nx ~ iiife.xPx. 



(13) 



We thus estimated the total number iV t e * t p of expected 
[IMBH,MSP] systems (i.e. those [IMBH,MSP] in which the 
radio beams of the MSP sweep the direction to the Earth), 
summing over all channels and over the sample of GCs host- 
ing at least one MSP. We find A t c xp ~ 0.1, if a ~ 100 M 
IMBH is hosted in all the GCs which are currently known 
to include a MSP. Thus, the detection of an [IMBH,MSP] 
binary has at present a low probability of occurrence 5 . 

The derived value of iV t c * t p is a firm lower limit since 
n-MSP represents a lower limit to the MSP density in a GC 
core, given that we considered only the already detected 
MSPs. The ongoing deep surveys running at GBT (Ransom 
et al. 2005), GMRT (Freire et al. 2004) and Parkes (Possenti 
et al. 2003) are rapidly increasing the known population of 
MSPs in GCs, suggesting that additional clusters may con- 
tain a rich population of MSPs. The likelihood of unveiling a 
binary [IMBH, MSP] will become significantly higher when 



6 DETECTABILITY OF [IMBH,MSP] 

BINARIES IN GLOBULAR CLUSTERS 

In Section 3 and 5 we investigated the formation of binaries 
hosting an IMBH and a MSP, via single-binary and binary- 
binary interactions. Here, we compute their formation rates 



5 No strong bias against the detection of an [IMBH, MSP] binary 
is caused by its the orbital motion. In fact, Patruno et al. (2005) 
showed that the discovery of a bright MSPs orbiting an IMBHs 
at mean separations of a few AU is not hampered by the Doppler 
modulation of the radio pulses. 
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new more powerful radio telescopes will become available. 
In particular the planned SKA (Cordes et al. 2004) is ex- 
pected to improve of 1-2 orders of magnitude the sensitivity 
limits of the present instruments. That will allow to probe 
the faintest end of the luminosity function of the MSPs in 
GCs. If the current extrapolations of this luminosity func- 
tion (Ransom et al. 2005; Camilo & Rasio 2005) will turn 
out to be correct, an order of magnitude more MSPs could 
be found in the core of the Galactic GCs, that have been 
missed by the current surveys due to their relative faint- 
ness. In this case, A^ t p « 1 and SKA will be able to detect 
all of this kind of systems. Therefore, a complete search for 
MSPs in the GCs of the Milky Way with SKA will have the 
potentiality of testing the hypothesis that IMBHs of order 
100 Mq are commonly hosted in GCs. 

The detection of one [IMBH,MSP] system will immedi- 
ately give the chance of measuring the mass of the IMBH 
from pulsar timing with at least 1% accuracy (Cordes et al. 
2004). Even more interesting, the presence of a very stable 
clock (like MSPs usually are) orbiting a probably rotating 
~ 100 M Q black hole makes this system a potentially unique 
laboratory of relativistic physics. In fact, many still elusive 
higher order relativistic effects depend on the spin and on 
the quadrupole moment of the rotating black hole (Wex & 
Kopeikin 1999) and the latter two quantities scale with the 
mass squared and the mass cubed of the BH, respectively. 
Therefore, an [IMBH, MSP] binary is a more promising tar- 
get for studying the physics in the surroundings of a BH 
(Kramer et al. 2004) than a binary comprising a MSP and 
a stellar mass BH. 



a MSP can allow for the formation of a stable hierarchical 
triple, where the MSP occupies the external orbit. When the 
internal [IMBH,BH] binary merges due to orbital decay by 
gravitational waves emission, the triple evolves into a new 
[IMBH,MSP] binary. 

• The [IMBH, MSP] binaries are expected to form with 
very high eccentricities (e ~ 0.9) and tight orbits (< 7 AU). 
This means that they can be important sources of gravi- 
tational waves, either in the in-spiral phase or in the final 
merging event. 

• Due to the aforementioned gravitational quadrupole ra- 
diation, the [IMBH, MSP] binaries are relatively short-lived, 
in-spiraling to coalescence in ~ 10 8 yr. This lifetime is sig- 
nificantly shorter than the estimated formation timescale 
of [IMBH, MSP] binaries which may be detectable with the 
present instrumentation. 

• If IMBHs of ~ 100 Mq are commonly hosted in the 
Galactic GCs, next-generation radio telescopes, like SKA, 
will have the possibility of detecting at least one of these 
exotic binaries. 
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7 SUMMARY 

In this paper, we investigated the dynamical processes lead- 
ing to the capture of a MSP by an IMBH in the dense core 
of a GC. We simulated single-binary and binary-binary en- 
counters between an IMBH and a MSP, either single or with 
a WD companion. The binary MSPs have masses and orbital 
parameters chosen according to the distribution observed in 
a sample of 23 GCs. In order to account for all the possible 
configurations of IMBHs hosted in GCs, we have consid- 
ered the case of a single IMBH, of an [IMBH, star] binary 
and of an [IMBH,BH] binary. For each of these cases we 
derived the cross-section for the formation of [IMBH, MSP] 
and [IMBH, WD] binaries, as well as the distribution of the 
final semi-major axes and eccentricities of such newly formed 
binaries. 

The main outcomes from this study are: 

• Dynamical encounters of a MSP with either single 
IMBHs or [IMBH, star] binaries promote the formation of 
[IMBH,MSP] binaries in ~ 10% and ~ 1 - 5% of the calcu- 
lated interactions, respectively. Similar rates were found for 
the formation of [IMBH,WD] binaries. The final distribu- 
tions of semi-major axes and eccentricities of the formed 
[IMBH,MSP] and [IMBH,WD] binaries are found to be 
in agreement with previous semi-analytical models (Pfahl 
2005). 

• We found that the presence of a stellar mass BH, or- 
biting around the IMBH, strongly inhibits the formation of 
an [IMBH, MSP] binary. Only in a small minority of cases 
(~ 0.2%), interactions between an [IMBH,BH] binary and 
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APPENDIX A: INITIAL CONDITIONS 

Al Initial semi-major axis distribution 

We describe here in some detail how we generate the initial 
distribution for the semi-major axis of the IMBH binaries. 

• [IMBH, star]: We have followed the analysis of Pfahl 
(2005) who considers the tidal disruption of a stellar binary 
off an IMBH. From considerations on energy conservation, 
the semi-major axes a* of the newly formed binary follows 
the relation 



1 m b /Mimbh\ 2/3 
a* = -a b I 

2 m CS c V mb / 



(Al) 



where m,b is the mass of the initial binary, ab its semi-major 
axis, and m csc the mass of the escaping star (see also the 
discussion in Section 3.2). To reproduce the initial distribu- 
tion for a*, we have considered a uniform distribution for 
the mass ratio of the stellar binary q = mi/m.2 and a distri- 
bution homogeneous in log(ab) for the values of the semi- 
major axes of the incoming binary in the range [0.01,10] 
AU. The upper and lower limits obtained are reported in 
Table 1. Fig. Al shows the initial distributions of a* (left 
panel) and e„ (right panel). Note that the distribution of 
a* is harder than that found in Blecha et al. (2006) 6 . If the 
real distribution would be less hard as in Blecha et al., our 
resulting [IMBH,MSP] formation rates for the [IMBH,star] 
case should be considered as a lower limit. Indeed a less 
bound initial companion to the IMBH would be more easily 
ejected by the unstable triple interaction with the MSP (see 
Section 3.2). 

• [IMBH,BH]: We expect that [IMBH,BH] binaries can form 
dynamically in the core of a GC. If the IMBH has been 
formed through a succession of gravitational encounters with 
stellar-mass BHs, then we expect some of these to be ejected 
in the outer region of the GC and to sink back to the core 
by dynamical friction (Sigurdsson & Hernquist 1993). The 
formation of the [IMBH,BH] binary can then be the result 
of one of the following interactions: 
IMBH+[BH,star] [IMBH,BH]+star; 
[IMBH,star] + [BH,star] -> [IMBH,BH]+stars; 
[IMBH,star] + BH [IMBH,BH]+star. 

The [IMBH,BH] binary just formed is assumed to have 
a separation comparable to the IMBH influence radius. This 
is not our initial condition for simulating the encounters 
with the [MSP, WD] systems, since we have accounted for 
the intrinsic long term evolution of the [IMBH,BH] binary 
parameters. Accordingly, we have generated the values of the 
initial [IMBH,BH] binary semi-major axis, (i.e. the values of 
ttBH from which we start the 3 or 4 body simulations) from a 
distribution obtained sampling uniformly in time when con- 
sidering the evolution of obh due to the hardening (i) off 
cluster stars, and (ii) by gravitational wave back-reaction. 

In phase (i), denoted in the text as [IMBH,BH]h,*, the 
evolution of obh is governed by the equation: 



6 We note that in their simulation, Blecha et al. consider stel- 
lar cluster considerably different from ours. Indeed, they study 
the formation of [IMBH, star] binaries in young clusters (their 
simulation stops after 100 Myr) with a correspondingly different 
population of stellar binaries. We argue that this can be the main 
cause of the difference in the distribution of a*. 
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Figure Al. Initial distribution of the semi- major axes (left 
panel) and eccentricities (right panel) of the initial states 
[IMBH.star] for the M IM bh = 100 M Q . 
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holding until a B H = a g w(eBH = 0.7) set by equation (3) 
(Hills 1975). In equation (A2) we assumed fixed the values 
of (p») = 7 x 10 5 M Q pc~ 3 and (v*) = lOkms -1 inferred av- 
eraging over the current GC sample described in Section 6. 
A change in and (t>») due to the internal evolution of 
the GC should also change the obh distribution. In particu- 
lar, a lower value for the stellar density should enhance the 
right end tail of the distribution. We argue that in this case 
the formation of [IMBH,MSP] binaries could be enhanced. 
Indeed, the presence of an initial companion, bound to the 
IMBH on a less tight orbit than that considered in our study, 
would be more easily ejected by the MSP (see Section 3.2). 

In phase (ii), the binary hardens by gravitational waves 
back-reaction (phase denoted with [IMBH,BH] gw ). The evo- 
lution of the orbital parameters are given by (Peters 1964): 

da B H 64 G 3 m B HMlMBH ( 

/ (6bh ) (A3) 
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Figure A2. Upper panel: the distribution for the initial ec- 
centricity of the [IMBH,BH] binaries are shown in the regimes 
of hardening off stars (right) and gravitational waves (left), for 
Mjmbh = 100 Mq. Lower panel: the distribution of the initial 
semi-major axes for the same systems are shown: the solid ver- 
tical line separates the gravitational wave regime (left) and the 
hardening off stars regime (right). 



de B H _ 304 G 3 m B HMiMBH (ra B H + Mimbh) 
d£ 15 c 5 ag H 

where 

/(e B H) = (l - e| H ) 7/2 ( 



73 2 37 4 
' + 24 6BH + 96 6BH 



Sr(e BH ) (A4) 



(A5) 



,-,/-. 2 \-5/2 /, 121 o \ 

9(eBH) = (1 - e BH ) e BH (1 + — e BH J 



(A6) 



dt 



c J a BH 



The above equations (A3-A6) are integrated with the initial 
condition: (cbh) and a trial distribution for eBH 

that follows the thermal distribution. Fig. A2 shows the re- 
sulting distribution for obh and eBH during the two different 
regimes. 7 



A2 The integration 

In this subsection we describe details on the integration of 
three-body and four-body encounters with the codes Chain 
and FEBO. For each run we divide the integration into two 
parts. (I) We consider the two binaries as point-like objects 
until their centers of mass are at a distance larger than 50 
times the semi-major axis of the IMBH binary. (II) When 
this critical distance is reached, we start the four- body 
integration. As a consequence, the time spent in the two 
body approximation decreases as the semi-major axes of the 



7 Note that in phase (ii), the distribution should not be affected 
by any change in the structural parameters of the GC, depending 
only on the orbital parameters (see equations (A3) and (A4)). 
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IMBH binary become wider. Correspondingly, the overall in- 
tegration time gets longer the wider the semi-major axis of 
the IMBH binary is, and it becomes prohibitively long for 
large values of obh (or a* for the [IMBH, star] binaries). For 
this reason, we insert a cut-off at 5 AU. For wider systems, 
we expect that the available binding energy of the IMBH bi- 
nary is insufficient to unbind the [MSP, WD] binary, so that 
the ionization of the binary can be mainly due to the tidal 
effect of the massive IMBH. 



A3 Impact parameters 

In this subsection we focus attention on the choice of the 
maximum impact parameter 6 max for a correct determina- 
tion of the cross section. According to gravitational focusing, 
a point mass with impact parameter b, moving at infinity 
with relative velocity Voo, has pericenter 



P ■ 



2GMn 



(A7) 



For the single IMBH, the maximum value of the peri- 
astron p m ax is set at a few tidal radii tt ; while for a binary 
IMBH, the value of the maximum impact parameter for a 
non negligible energy exchange is typically limited up to a 
value of the order of a few semi-major axis of the binary 
IMBH, i.e. p max ~ xcibh or p m ax ~ xa* (Hills 1983), where 
x is close to 3 in all cases. In each run 6^ is assigned using 
equation (A7). In order to guarantee that we have accounted 
for all the impact parameters leading to the formation of 
an [IMBH,MSP] (or [IMBH, WD]) binary, and to guarantee 
cross section convergence, we verified a posteriori that the 
distribution of all relevant b leading to the desired end- 
states, drops to zero well before i>max- Fig. A3 illustrates the 
case of the single 100 M IMBH interacting with the PSR- 
A like [MSP, WD] binary. The encounters ending with the 
formation of [IMBH, MSP] binaries are the ones represented 
in the hatched area. Clearly, the distribution drops to zero 
well before &max- 

In the channels where either the tidal radius of the in- 
coming binary MSP, or the semi-major axis of the IMBH 
binary vary from encounter to encounter (according to the 
initial distributions described in Sections 2.2 for the binary 
IMBH, and in Section 5 for the binary MSP), we allowed 
frmax t° var y accordingly, and defined a mean (6^ ax ), ob- 
tained averaging over all choices of rr and/or a* (abh). This 
average is used to compute the cross section in equation (5). 
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Figure A3. Distribution of impact parameters giving rise to the 
formation of a [IMBH, MSP] system (hatched histogram), com- 
pared to the initial (empty histogram) for the case of encounters of 
PSR-A like [MSP, WD] binaries off the single 100 M IMBH. The 
hatched distribution drops to zero at 300 AU 2 , while 6^ = 1000 
AU 2 . 



